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Abstract 

Wightman function, the vacuum expectation values of the field square and the energy- 
momentum tensor are investigated for a massive scalar field with general curvature coupling 
parameter in the region between two coaxial cylindrical boundaries. It is assumed that the 
field obeys general Robin boundary conditions on bounding surfaces. The application of a 
variant of the generalized Abel-Plana formula allows to extract from the expectation values 
the contribution from single shells and to present the interference part in terms of exponen- 
tially convergent integrals. The vacuum forces acting on the boundaries are presented as the 
sum of self-action and interaction terms. The first one contains well-known surface diver- 
gences and needs a further renormalization. The interaction forces between the cylindrical 
boundaries are finite and are attractive for special cases of Dirichlet and Neumann scalars. 
For the general Robin case the interaction forces can be both attractive or repulsive depend- 
ing on the coefficients in the boundary conditions. The total Casimir energy is evaluated by 
using the zeta function regularization technique. It is shown that it contains a part which is 
located on bounding surfaces. The formula for the interference part of the surface energy is 
derived and the energy balance is discussed. 



PACS numbers: ll.lO.Kk, 03.70. -Fk 



1 Introduction 

The Casimir effect is one of the most interesting macroscopic manifestations of the nontrivial 
structure of the vacuum state in quantum field theory. The effect is a phenomenon common 
to all systems characterized by fluctuating quantities and results from changes in the vacuum 
fluctuations of a quantum field that occur because of the imposition of boundary conditions or 
the choice of topology. It may have important implications on all scales, from cosmological to 
sub nuclear, and has become in recent decades an increasingly popular topic in quantum field 
theory. In addition to its fundamental interest the Casimir effect also plays an important role 
in the fabrication and operation of nano- and micro-scale mechanical systems. The imposition 

* E- mail : sahary an® server . phy sdep . r . am 



1 



of boundary conditions on a quantum field leads to the modification of the spectrum for the 
zero-point fluctuations and results in the shift in the vacuum expectation values for physical 
quantities such as the energy density and stresses. In particular, the confinement of quantum 
fluctuations causes forces that act on constraining boundaries. The particular features of the 
resulting vacuum forces depend on the nature of the quantum field, the type of spacetime 
manifold, the boundary geometries and the specific boundary conditions imposed on the field. 
Since the original work by Casimir U many theoretical and experimental works have been 
done on this problem (see, e.g., [HISlllllS] and references therein). Many different approaches 
have been used: mode summation method with combination of the zeta function regularization 
technique. Green function formalism, multiple scattering expansions, heat-kernel series, etc. 
Advanced field-theoretical methods have been developed for Casimir calculations during the 
past years jO]. 

An interesting property of the Casimir effect has always been the geometry dependence. 
Straightforward computations of geometry dependencies are conceptually complicated, since 
relevant information is subtly encoded in the fluctuations spectrum. Analytic solutions can usu- 
ally be found only for highly symmetric geometries including planar, spherically and cylindrically 
symmetric boundaries. Aside from their own theoretical and experimental interest, the problems 
with this type of boundaries are useful for testing the validity of various approximations used 
to deal with more complicated geometries. In particular, a great deal of attention received the 
investigations of quantum effects for cylindrical boundaries. In addition to traditional problems 
of quantum electrodynamics under the presence of material boundaries, the Casimir effect for 
cylindrical geometries can also be important to the flux tube models of confinement [ZllHl and for 
determining the structure of the vacuum state in interacting field theories . The calculation 
of the vacuum energy of electromagnetic field with boundary conditions defined on a cylinder 
turned out to be technically a more involved problem than the analogous one for a sphere. First 
the Casimir energy of an infinite perfectly conducting cylindrical shell has been calculated in 
Ref. ^01 by introducing ultraviolet cutoff and later the corresponding result was derived by zeta 
function technique jlll 1121 [T^ . The local characteristics of the corresponding electromagnetic 
vacuum such as energy density and vacuum stresses are considered in jjl] for the interior and 
exterior regions of a conducting cylindrical shell, and in jl5| for the region between two coaxial 
shells (see also |16j). The vacuum forces acting on the boundaries in the geometry of two cylin- 
ders are also considered in Refs. jl7j . Less symmetric configurations of a semi-circular infinite 
cylinder and of a wedge with a coaxial cylindrical boundary are investigated in Refs. |18| I19j. 
A large number of papers is devoted to the investigation of the various aspects of the Casimir 
effect for a dielectric cylinder (see, for instance, and references therein). From another per- 
spective, the influence of a dielectric cylinder on the radiation process from a charged particle 
has been discussed in pT] . 

An interesting topic in the investigations of the Casimir effect is the dependence of the 
vacuum characteristics on the nature of boundary conditions imposed. In Ref. 1^ scalar 
vacuum densities and the zero-point energy for general Robin condition on a cylindrical surface 
in arbitrary number of spacetime dimensions are studied for massive scalar field with general 
curvature coupling parameter. The Robin boundary conditions are an extension of the ones 
imposed on perfectly conducting boundaries and may, in some geometries, be useful for depicting 
the finite penetration of the field into the boundary with the 'skin-depth' parameter related to 
the Robin coefficient [2'A\ I24j . It is interesting to note that the quantum scalar field satisfying 
Robin condition on the boundary of a cavity violates the Bekenstein's entropy-to-energy bound 
near certain points in the space of the parameter defining the boundary condition [2S1- This 
type of conditions also appear in considerations of the vacuum effects for a confined charged 
scalar field in external fields [201 ™ quantum gravity [23 12B1 120] • Mixed boundary conditions 
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naturally arise for scalar and fermion bulk fields in braneworld models [30] . For scalar fields with 
Robin boundary conditions, in Ref. '31] it has been shown that in the discussion of the relation 
between the mode-sum energy, evaluated as the sum of the zero-point energies for each normal 
mode of frequency, and the volume integral of the renormalized energy density for the Robin 
parallel plates geometry it is necessary to include in the energy a surface term concentrated on 
the boundary (for the discussion of the relation between the local and global characteristics of 
the vacuum see also Ref. |32( I33| I34| \22\ I35j ) . An expression for the surface energy- momentum 
tensor for a scalar field with a general curvature coupling parameter in the general case of bulk 
and boundary geometries is derived in Ref. ."}() . Related cosmological constant induced on the 
brane by quantum fluctuations of a bulk field in braneworld scenarios has been considered in 



In the present paper, we consider the Casimir densities in the region between two coaxial 
cylindrical shells on background of the {D + l)-dimensional Minkowski spacetime. The posi- 
tive frequency Wightman function, the vacuum expectation values of the field square and the 
energy-momentum tensor are investigated for a massive scalar field with general curvature cou- 
pling parameter. In addition to describing the physical structure of the quantum field at a given 
point, the energy-momentum tensor acts as the source of gravity in the Einstein equations. It 
therefore plays an important role in modelling a self-consistent dynamics involving the gravita- 
tional field P3 . For the general case of Robin boundary conditions with different coefficients for 
the inner and outer boundaries, we derive formulae for the forces acting on the boundaries due 
to the modification of the spectrum of the zero-point fiuctuations by the presence of the second 
boundary. The Casimir energy and the surface energy are investigated as well and the energy 
balance is discussed. 

The plan of the paper is as follows. In the next section we derive a formula for the Wight- 
man function in the region between two cylindrical surfaces. The reason for our choice of the 
Wightman function is that this function also determines the response of the particle detectors in 
a given state of motion. To evaluate the bilinear field products we use the mode-sum method in 
combination with the summation formula from I34j . These formula allows (i) to extract from 
vacuum expectation values the parts due to a single cylindrical boundary, and (ii) to present the 
interference parts in terms of exponentially convergent integrals involving the modified Bessel 
functions. The vacuum expectation values of the field square and the energy-momentum tensor 
are obtained from the Wightman function and are investigated in section |3J The vacuum forces 
acting on the bounding surfaces are considered in section HI They are presented as the sum of 
self-action and interaction terms. The formulae are derived for the interaction forces between 
the cylinders. Section |51 is devoted to the total vacuum energy evaluated as a sum of zero-point 
energies of elementary oscillators. We show that this energy in addition to the volume part 
contains a part located on the bounding surfaces. The formula for the interaction part of the 
surface energy is derived. Further we discuss the relation between the vacuum energies and 
forces acting on the boundaries. Section |21 concludes the main results of the paper. 

2 Wightman function 

Consider a real scalar field ip with curvature coupling parameter ^ satisfying the field equation 



where R is the curvature scalar for a (D + l)-dimensional background spacetime, Vj is the 
covariant derivative operator. For special cases of minimally and conformally coupled scalars 
one has ^ = and = = (D — 1)/4Z), respectively. Our main interest in this paper will be 
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one-loop quantum vacuum effects induced by two infinitely long coaxial cylindrical surfaces with 
radii a and b, a < b, in the Minkowski spacetime. For this problem the background spacetime is 
flat and in Eq. we have R = 0. As a result the eigenmodes are independent on the curvature 
coupling parameter. However, the local characteristics of the vacuum such as energy density 
and vacuum stresses depend on this parameter. In accordance with the problem symmetry we 
will use cylindrical coordinates (r, (j),zi, zn), N = D — 2, and will assume that the field obeys 
Robin boundary conditions on bounding surfaces: 

(^Aj + S,nJ^)Vi) if (x) l^. = 0, i = a, b, (2) 

with Aj and Bj being constants, re^^-^ is the inward-pointing normal to the bounding surface 

r = j. For the region between the surfaces, a ^ r ^ b, one has n^^.^ = nj5\ with the notations 
Ua = 1 and Hf, = —1. The imposition of boundary conditions on the quantum field modifies the 
spectrum for zero-point fluctuations and leads to the modification of the vacuum expectation 
values (VEVs) for physical quantities compared with the case without boundaries. First we 
consider the positive frequency Wightman function. The VEVs of the field square and the 
energy-momentum tensor can be evaluated on the base of this function. By the same method 
described below any other two-point function can be evaluated. 

Let (x)} is a complete orthonormal set of positive and negative frequency solutions 

to the field equation, specified by a set of quantum numbers a and satisfying the boundary 
conditions By expanding the field operator and using the standard commutation relations, 
the positive frequency Wightman function is presented as the mode-sum 

{0\^{xMx')\0) = Y,Va{x)^*a{x), (3) 

a 

where |0) is the amplitude for the corresponding vacuum state. For the region a ^ r ^ b the 
eigenfunctions are specified by the set of quantum numbers a = (n, 7, k), n = 0, ±1, ±2, . . ., and 
have the form 

(Pa{x) = I3ag\n\{'ia,-fr) exp{in(j) + ikrii - iujt), (4) 

with ry = (zi, Z]\f), LO = v/+^^) = + '^^j s-'^d 

gnha,^r) = Y^^^^apnM " Ji''\ia)Ynhr). (5) 

In Eq. Jniz) and Yn{z) are the Bessel and Neumann functions and for a given function f{z) 
we use the notation 

P\z) EE A,f{z) + {njB,/j)zf'{z), (6) 

with j = a, b. The eigenfunctions (j!]) with the radial part © satisfy the boundary condition 
on the inner surface. The eigenvalues for the quantum number 7 are quantized by the boundary 
condition Q on the surface r = b. Prom this condition it follows that the possible values of 7 
are solutions to the equation (r/ = b/a) 

C^\v,ia) ^ jt\7a)Yj^'\^b) - Yt\^a)J^\^b) = 0. (7) 

In the discussion below the corresponding positive roots we will denote by 7a = On^u / = 1, 2, . . ., 
assuming that they are arranged in the ascending order, an,i < Cn.i+i- 

The coefficient Pa in Q is determined from the orthonormality condition for the eigenfunctions 

J dV(paix)ipl,{x) = -^dnn'^ll'^O^ - k'), (8) 
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where the integration goes over the region between the cyhndrical shells. By making use of the 
standard integral for the cylindrical functions (see, for instance, |39j^. one finds 



Pi 



4a;a(27r)^-i ' 



(9) 



with the notation 



-1 



(10) 



Substituting eigenfunctions ^ into the mode-sum formula Q, for the positive frequency 
Wightman function one finds 



oo oo 



(oi^(xMx')io> = ^ d^k^'Y: 



zgn{z, zr/a)gn{z, / ^^Ij^ab ^ ^ 



X cos [n (</) — 0')] exp[— ( ry — rji 



{t-t')] 



where the prime on the summation sign means that the summand with n = should be halved. 
As the expressions for the eigenmodes an,i are not explicitly known, the form for the 
Wightman function is inconvenient. For the further evaluation of this VEV we apply to the sum 
over / the summation formula |161 134j 



1=1 



oo 

I 



h{x)dx 



{a)2. 



IT 



dxnan{x,rix) [/i(xe^*/2) + h{xe-^'/^] 



(12) 



where 



K^\y)/Kr{x 



(13) 



and Inix), Kn{x) are the Bessel modified functions. Here we have assumed that all zeros for 
the function C^{rj, z) are real. In case of existence of purely imaginary zeros we have to include 
additional residue terms on the left of formula p2|) (see ^Mi ) • Formula p2j) is valid for functions 
h{z) satisfying the condition \h{z)\ < e(x)e^l^l, z = x + iy, c < 2{rj — 1), for large values \z\, 
where x'^^^°-^~^e{x) for x oo, and the condition h{z) = o{z^^) for z — > 0. 
For the evaluation of the Wightman function, as a function h[x) we choose 



,/ \ xgn{x,xr/a)gn{x,xr /a) / o / i I PTw j.m 



(14) 



The corresponding conditions are satisfied if r + r' + \t — t'\ < 2b. In particular, this is the 
case in the coincidence limit t = t' for the region under consideration. Now we can see that the 
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application of formula allows to present the Wightman function in the form 



[0\ip{x)ip{x')\0) 



(2vr) 



^^'cos [n 



71=0 



*')]/ 



2 /■ xfiai^faz, 6z) 



xG(f)(az, zr)G(,")(az, zr') cosh {t - t')y^z^-kl 

with the notations (the notation with j = b will be used below) 

G(f)(x,y) = Ki^\x)Uy) - li^\x)Kn{y), j = a,b. 



(15) 



(16) 



In the limit b ^ oo the second term in figure braces on the right of (|15p vanishes, whereas 
the first term does not depend on b. It follows from here that the part with the first term 
presents the Wightman function in the region outside of a single cylindrical shell with radius a 
(of course, this may also be seen by direct evaluation of the corresponding Wightman function). 
To simplify this part we use the identity 



gniz,zr/a)gn{z,zr' /a) 



Jt^\z)+Yr'iz) 



1 ^ j''"'\z) 
Jn{zr/a)Jn{zr'/a) - - ^ 

xHt\zr/a)Ht\zr'/a), 



(17) 



with H^\z), a = 1,2 being the Hankel functions. Substituting this into the first integral in the 
figure braces of Eq. H15|l we rotate the integration contour over z by the angle 7r/2 for a = 1 and 
by the angle —tt/2 for a = 2. Under the condition r + r' — |t — t'| > 2a, the integrals over the 
arcs of the circle with large radius vanish. The integrals over {0,iakm) and {0,—iakm) cancel 
out and after introducing the Bessel modified functions one obtains 



dz 



, = /^,/n(^r)J„(zrO ^ 

J<")^(z) + Fi'^)^(z) / v/i^Tfcl 

- ^ /d.z4T^^^^^?^^^-«My^^(i'-o].(i8) 



Substituting this into formula ()15() . the Wightman function is presented in the form 

{0\^ixMx')\0) = {ipixMx')f^ + {^{xMx'))^"^ 

22-D 

n=0 



CO ^ 

?^^'cos [n{cp-^')] 1 AeK''"-'"'l) J 



dz z 



Qgnjaz, bz) 
\J z^ — k"^ 



xG^^\az, zr)G^^\az, zr') cosh \{t - t')yG^^M 



(19) 



where {<p{x)ip(x'))^^'^ is the Wightman function for a scalar field in the unbounded Minkowskian 
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spacetime, and 



oo 



^ n=0 



fn\az) Kn{zr)Kn{zr') 



X / dz z -\ -ny -n^' > cosh[v^i^3^(t' - t)] (20) 

is the part of the Wightman function induced by a single cylindrical shell with radius a in the 
region r > a. Hence, the last term on the right of (|19|) is induced by the presence of the second 
shell with radius b. 

By using the identity 

njQ,jn{az, hz)G^^\az, zr)G[^\az, zr') 

j=a,b 

= ^$^In{zr)Uzr') - ^^^K,Xzr)Kr.{zr'l (21) 
In [bz) K\'{az) 

with the notation 

the Wightman function can also be presented in the equivalent form 

{0\^ixMx')\0) = {ip{xMx')f^ + {^ixMx')f^ 

oo 



71=0 



xG'^\bz,zr)G'^\bz,zr') cosh {t - t')yj z^ - kl^ , (23) 

where 



{^{xMx')f^ = -^X;'cos[n(<A-</<')] /rf^ke^^ 



n=0 

oo 



X j dzz ^l}^''^ In{ zr)Uzr ') ^ _ 

/Cm 

is the part induced by a single cylindrical shell with radius b in the region r < b. Note that 
formulae (|2U|1 and (|24j) are related by the interchange a b, In ^ Kn- In the formulae for the 
Wightman function given above the integration over the angular part of the vector k can be 
done with the help of the formula 

e^'^-Fik) r •^^/2-i(^|x|) 



d-k—-^= dkk^'-'F{k) ;: 7:^/ ,^ ' (25) 

(27r)- JQ (A;|x|)^/^ ^ 



for a given function F{k). 
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3 VEVs of the field square and the energy-momentum tensor 



Having the Wightman function we can evaluate the VEVs of the field square and the energy- 
momentum tensor. These VEVs in the regions r < a and r > b are the same as those for a 
single cylindrical surface with radius a and b respectively and are investigated in [^. For this 
reason in the discussion below we will be concerned on the region a < r < b. By making use 
of formulae H19|) and H23|) for the Wightman function and taking the coincidence limit of the 
arguments, for the VEV of the field square one finds 



2\(0) 



(Olv^'lO) = {if 



where j = a and j 



2\{i) 



n=0 



AdY^ I duu{u^ -rr?)~ ^jn{auM)G'i^'^{ju,ru), (26) 



b provide two equivalent representations and 

22-D 



A 



D 



TT 2 



T(^) 



(27) 



To obtain this result we have employed the integration formula 



dk I dz 



f{z) V^r{§-i 



1,2 



2r 



du (i 



m' 



(28) 



For points away from the boundaries the last two terms on the right of formula H26() are finite 
and, hence, the subtraction of the Minkowskian part without boundaries is sufficient to obtain 
the renormalized value for the VEV: (y?^)ren = (Olv^^lO) — {lp^)'^^\ In formula (|26|) the part (<^^)''-'^ 
is induced by a single cylindrical surface with radius j when the second surface is absent. The 
formulae for these terms are obtained from H20|) and 1)241) in the coincidence limit. For j = a 
one has 



2\(a) 



-A 



n=0 



du U 



U 



-D-3 
2 



■Kliru) 



(29) 



and the formula for (v?^)^^-* is obtained from here by the replacements a ^ b, I ^ K . The last 
term on the right of formula (|26|) is induced by the presence of the second cylindrical surface. 
The surface divergences in the boundary induced parts in the VEV of the field square are the 
same as those for single surfaces and are investigated in j^. In particular, the last term on the 
right of ()26j) is finite for r = j. It follows from here that if we present the renormalized VEV of 
the field square in the form 



2\(af.) 



(30) 



then the interference part {(p"^)^""^^ is finite on both boundaries. In the limit a ^ this part 
vanishes as a for Robin boundary condition on the inner shell {Aa,Ba 7^ 0), as for Neumann 
boundary condition (Aa = 0), and like 1/lna for Dirichlet boundary condition {Ba = 0). In 
the limit b ^ 00 and for a massless field the interference part vanishes as In 6/6^^^ for Robin 
boundary condition on the inner shell, as \nb/b^~^^ for Neumann boundary condition, and like 
l/b^~^ for Dirichlet boundary condition. In the same limit under the condition mb S> 1 the 
interference part is exponentially suppressed. 

The VEV for the energy-momentum tensor is obtained by using the formulae for the Wight- 
man function and the VEV of the field square: 



{0\Tik\0) = lim did'k{0\^ix)ipix')\0) + 



(0|^». (31) 
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Substituting (|19|). (|26l) into l|31|) one finds (no summation over i) 



{0\Tt\0) = {Tt)^'H{Th^^^ + AnS'^Y! 

n=0 



where we have introduced notations 



Fj^'^fiz)] 



1 - rr?r^ jz^ 



D - 1 



1 



n 



n 



+ ^f{z)f'iz), 



(32) 

(33a) 
(33b) 



^ 2 



f'Hz) + ( ^ + 1 ) fiz 



n 



-^f{z)f'{z) + ^f{z), (33c) 



with f{z) = Gn {ju,z) and 



Fj:nfiz)]=Fl,'nf{z)], i = 3,...,D-l. 



(34) 



In formula (|32)) the term (Tj^)^^^ is induced by a single cylindrical surface with radius j. These 
parts for both interior and exterior regions are investigated in [22 . The formula for the case 
j = a is obtained from H29|) by the replacement K^{ru) — > u'^Fn\Kn{ru)]. As in the case of 
the field square, the renormalized VEV of the energy-momentum tensor can be presented in the 
form 

{Tt)ren=Y.(^')^'^ + (^')^^'^' (35) 
j=a,b 



where the surface divergences are contained in the single boundary parts only and interference 
part is finite on the boundaries. The explicit formula for the latter is obtained by subtracting 
from the last term on the right 1)32(1 the corresponding single surface part. Two equivalent 
representations are obtained by taking in ()32|) j = a or j = b. Due to the presence of boundaries 
the vacuum stresses in radial, azimuthal and axial directions are anisotropic. For the axial stress 
and the energy density we have standard relation for the unbounded vacuum. It can be easily 
checked that the separate terms in formulae (|26() and (|32() satisfy the standard trace relation 



T;* = ^(c-eD)v,vv' + mV, 



(36) 



and the continuity equation VjT|, 
the form 



0. For the geometry under consideration the latter takes 



At/ + i {Tl 

ar r 



r|) = 0. 



(37) 



In particular, this means that the r-dependence of the radial pressure leads to the anisotropy of 
the vacuum stresses. For a conformally coupled massless scalar the vacuum energy- momentum 
tensor is traceless. In the limit a, 6 — > c« with fixed b — a from the formulae above the results for 
the geometry of two parallel plates with Robin boundary conditions are obtained. In the limit 
a — ^ and for Aa ^ the main contribution into the interference part of the vacuum energy- 
momentum tensor comes from n = term and this part behaves like a in the case Ba ^ and 
like 1/ ln(a/6) in the case Ba = 0. In the same limit and for Aa = the interference part behaves 
as o^. For a massless scalar field in the limit 6 — > oo and for Aa, Ba ^ the main contribution 
comes, again, from n = term. The corresponding energy density behaves as b^~^ , while the 



9 



vacuum stresses behave like In 6/6^^^. For a massive scalar field under the condition mb ^ 1 
the interference part of the energy-momentum tensor is exponentially suppressed. 

The interference part of the vacuum energy-momentum tensor gives finite contribution into 
the vacuum energy in the region a ^ r ^ 6. To evaluate this contribution we note that the 
function Gn\ju, ru) is a modified cylindrical function with respect to the argument r. By using 
the formula for the integral involving the square of a modified cylindrical function (see, for 
instance, (39j) and noting that for these functions 



z 



the following formula can be obtained 

rb 



2^ / drr{T^,)^^') = eI%^^ + + Ai?(-^ (39) 

J a 

where E^^^^^ (-^o<r<a) vacuum energy in the region h ^ r < co (0 ^ r ^ a) for a single 

cylindrical shell with radius a (b) and 



/ D — 3 

A£;(^°^) = ttAdY^ / duu{u^ -m^)~ ^ njnjn{au,bu) 

n=0 -^"^ j=a,b 

X |(4^ - l)n,A,^ - ' {B]{u^ + n^f) - A^] } . (40) 



Now, by taking into account formula H39() . for the total volume energy in the region a ^ r ^ b 
one finds 



J a 

where E^'^"^^ (-^r <I°'^ ) volume part of the vacuum energy outside (inside) a single cylin- 

drical shell with radius a (6). Of course, due to the surface divergences these single shell parts 
cannot be obtained directly by the integration of the corresponding densities and need additional 
renormalization. As we will see below in section |3 the total vacuum energy in addition to the 
volume part contains the contribution located on the boundaries. 



4 Interaction forces 

Now we turn to the vacuum forces acting on the cylindrical surfaces. The vacuum force per 
unit surface of the cylinder at r = j is determined by the } - component of the vacuum energy- 
momentum tensor at this point. For the region between two surfaces the corresponding effective 
pressures can be presented as the sum of two terms: 

P^'^=P?+4ly j = a,b. (42) 

The first term on the right is the pressure for a single cylindrical surface at r = j when the second 
surface is absent. This term is divergent due to the surface divergences in the vacuum expectation 
values and needs further renormalization. This can be done, for example, by applying the 
generalized zeta function technique to the corresponding mode-sum. This procedure is similar 
to that used in Ref. P^l for the evaluation of the total Casimir energy for a single cylindrical shell. 
This calculation lies on the same line with the evaluation of the surface Casimir densities and 
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will be presented in the forthcoming paper [30]. Here we note that the corresponding quantities 
will depend on the renormalization scale and can be fixed by imposing suitable renormalization 
conditions which relates it to observables. The second term on the right of Eq. (|42|) is the 
pressure induced by the presence of the second cylinder, and can be termed as an interaction 
force. Unlike to the single shell parts, this term is free from renormalization ambiguities and is 
determined by the last term on the right of formula 1)32(1 . Substituting into this term r = j and 
using the relations 

Gli^ {u, u) = -njB,/j, G^jy (n, n) = Aj/u, (43) 
for the interaction parts of the vacuum forces per unit surface one finds 



^(int) 



/ / D—3 



X [(nVj' + u^) B] + 4.^njAjBj/j - A]] . 



n=0 



(44) 



The expression on the right of this formula is finite for all non-zero distances between the 
shells. By taking into account the inequalities Kn{u)In{v) — Kn{v)In{u) > and K'^{u)I'^{v) — 
K'^{v)I'^{u) < for n < u, it can be seen that the vacuum effective pressures are negative 
for both Dirichlet and Neumann scalars and, hence, the corresponding interaction forces are 
attractive. For the general Robin case the interaction force can be either attractive or repulsive in 
dependence on the coefficients in the boundary conditions (see the numerical example presented 
below in figure [2|). The quantity determines the force by which the scalar vacuum acts on 
the cylindrical shell due to the modification of the spectrum for the zero-point fluctuations by 
the presence of the second cylinder. As the vacuum properties depend on the radial coordinate, 
there is no a priori reason for the interaction terms (and also for the total pressures p^^^) to be 
the same for j = a and j = b, and the corresponding forces in general are different. Note that 
the interaction parts act on the surfaces r = a+0 and r = 6 — 0. The vacuum forces acting on the 
sides r = a — and r = b + are the same as those for single surfaces. In combination with the 
parts p^^^ from (|42l) . the latter give the total vacuum forces acting on a single cylindrical shell. 
For Dirichlet and Neumann boundary conditions these forces can be obtained by differentiation 
of the corresponding Casimir energy (see below section For D = 3 massless scalar these 
forces are repulsive for Dirichlet case and are attractive for Neumann case |12j . 

Using the Wronskian for the Bessel modified functions, it can be seen that for j = a,b one 

has 

d 

Bjn{u)VLjn{au, bu) = jrij- 



dj 



In 



1 



lt\au)Ki'\bu) 



A'\bu)Kr{au 



(45) 



where na = 1, = — 1 and we have introduced the notation 

Bjniu) = B]{nyf + u^) + njAjBj/j - A]. 

This allows us to write the expressions ()44|) for the interaction forces per unit surface in another 
equivalent form: 



(46) 



„0) 



2j 



E' 

n=0 



duuiu — m 



■i^ 2 



X 



i + (4e-i)^|^ 



d_ 

dj 



In 



li^\au)Ki'\bu) 



fn\bu)Kl^\au) 



(47) 



This form will be used in the next section in the discussion of the relation between the vacuum 
forces and the Casimir energy. Note that for Dirichlet and Neumann scalars the second term on 
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the square brackets vanish and the interaction forces do not depend on the curvature couphng 
parameter. 

Now we turn to the investigation of the vacuum interaction forces in various hmiting cases. 
First of all let us consider the case when the radii of cylindrical surfaces are close to each other: 
h/a — 1 <^ 1. Noting that in the limit b ^ a the vacuum forces diverge and in the limit 
under consideration the main contribution comes from large values n, we can use the uniform 
asymptotic expansions for the Bessel modified functions (see, for instance, ST]). As the next 
step we introduce a new integration variable x 



m 



2 and replace the summation over n 
by the integration: X^n^o ~^ /o°° Further introducing a new integration variable n = ay and 
passing to the polar coordinates in the plane (x, y) after the integration of the angular part one 
finds 



\mt) 



vr#r(#) 



m 



D_ 
I 2 " 



{Aa-Bat){Ak-Bi,t) ^ 
{Aa + Bat){Ai, + Bkt)' 



-1 



(48) 



The latter formula coincides with that for the interaction forces between two parallel plates with 
Robin boundary conditions © on them j^. In this case the interaction forces are the same 
for both plates and do not depend on the curvature coupling parameter. Note that in the limit 
a ^ b with fixed values of the boundary coefficients and the shell radii, the renormalized single 
surface parts p^^^ remain finite while the interaction part goes to infinity. This means that for 
sufficiently small distances between the boundaries the interaction term on the right of formula 
()42() will dominate. 

For small values of the ratio a/6 we introduce in (|44|1 a new integration variable bu = x and 
expand the integrand by using the formulae for the Bessel modified functions for small values 
of the argument. For Aa,Ba ^ in the leading order the main contribution comes from n = 
term and we have 



(a) 
P{mt) 



P. 



{b) 
(int) 



aBa 
^b^Ba 



duu (u 



duu (u 



D-3 K 



(b) 



{bu) 



jz_3 ^2^2 _ 4^A,Bh/b - Al 



(49) 
(50) 



In the same limit and for Bg = the corresponding asymptotic formulae for the interaction 
forces are obtained by the replacement S,Aa/Ba — > — l/[4aln'^(a/6)] in formula H49|) and by the 
replacement Aga/Ba — > 1/ In (a/6) in formula (|50j) . In the case of Neumann boundary condition 
on the surface r = a {Aa = 0) the main contribution into the vacuum interaction forces comes 
from n = and n = 1 terms with the leading behavior 



P, 



(int) 



ib) 

P{mt) 



A 



A „2 /-oo 

Apa 



462 



du u (u 



m 



n=0 



D-3 . ,/ 



li'Hbu) 



n=0 



{u^ + nyb'')B^-AU,B,/b-A- 
li'^\bu) 



(51) 



(52) 



In figure ^ we have plotted the dependence of the interaction forces per unit surface on the 
ratio of the radii for the cylindrical shells, a/6, for the cases of D = 3 Dirichlet (left panel) and 
Neumann (right panel) massless scalars. As we have mentioned before, the interaction forces 
in these cases are attractive. In figure |21 we present the same graphs for the Robin boundary 
conditions with the coefficients Ba = 0, Bh/(Abb) = —0.2. The left panel corresponds to the 
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Figure 1: Vacuum interaction forces acting per unit surfaces of the cylindrical shells as functions 
of a/h. The left panel corresponds to the case of D = 3 massless Dirichlet scalar and the right 
one is for Neumann scalar. 

minimally coupled scalar and the right one is for a conformally coupled scalar. As we see, for the 
considered example the interaction force are repulsive for small distances between the surfaces 
and are attractive for large distances. This provides a possibility for the stabilization of the radii 
by the vacuum forces. However, it should be noted that to make reliable predictions regarding 
quantum stabilization, the renormalized single shell parts p\ also should be taken into account. 




alb alb 



Figure 2: The same as in figurenfor minimally (left panel) and conformally (right panel) coupled 
Robin scalars with Ba = 0, Bi,/{Aib) = -0.2. 

5 Casimir energy 

In this section we consider the total vacuum energy for the configuration of two coaxial cylindrical 
boundaries. In the region between the boundaries the total vacuum energy per unit hypersurface 
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in the axial direction is the sum of the zero-point energies of elementary oscillators: 



(53) 

n=0 1=1 

The expression on the right is divergent and to deal with this divergence we take the zeta 
function approach (for the application of the zeta function technique to the calculations of the 
Casimir energy see j42| and references therein). We consider the related zeta function 



C{s) = M 



s+l 



s/2 



(54) 



where the parameter /i with dimension of mass is introduced by dimensional reasons. Evaluating 
the integral over k we present this function in the form 



C(^) 



with the partial zeta function 



^^+1 r(^) 



(4vr)f r(f)a 



N-s 



Y^'Cnis-N), 



(55) 



n=0 



2„2^-s/2 



m a 



(56) 



1=1 



We need to perform the analytic continuation of the sum on the right of (|55|) to the neighborhood 
of s = —1. An immediate consequence of the Cauchy's formula for the residues of a complex 
function is the expression 

Cnis) = ^ I dz(z2 + mV)-^/2^lnCf(r?,z), (57) 

where C is a closed counterclockwise contour in the complex z plane enclosing all zeros an,i- We 
assume that this contour is made of a large semicircle (with radius tending to infinity) centered 
at the origin and placed to its right, plus a straight part overlapping the imaginary axis and 
avoiding the points ziziam by small semicircles in the left half-plane. When the radius of the 
large semicircle tends to infinity the corresponding contribution into C„ (s) vanishes for Re s > 1. 
Let us denote by and the upper and lower halves of the contour C. The integral on the 
right of Eq. H57|) can be presented in the form 



Cn {s) 



1 

27ri 



c 



oz 



+ 



1 

27ri 



E 

a=l,2 



z~-j^\nz) 
d 



dz (z2+^2^2)-s/2_i^ 

(ja OZ 



+ —y I dz(z2+mV)-^/2|-ln 
27r« ^-^2 Jc'- dz 



_.J^\z)H^%z)_ 
Sf'HzUi'Hrjz)} 



(58) 



where H^\z) are the Hankel functions. After parameterizing the integrals over imaginary axis 
we see that the parts of the integrals over (0, zizima) cancel and we arrive at the expression 



Cn (s) = - sin ^ / dz {z 



IT 



2 ^2^2)-./2|. 1 1^ 

oz 



+ \n. 



z-Ki^\z) 



lt\z)K^"'{vz) 



+ ln 1 



Kt\z)irinz) 



(59) 
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The integral with the last term in square brackets on the right of this formula is finite at 
s = —{N + 1) and vanishes in the limits a — > or 6 — > cxd. The integrals with the first and 
second terms in the square brackets correspond to the partial zeta functions for the region 
inside a cylindrical shell with radius b and for the region outside a cylindrical shell with radius 
a, respectively. As a result the total energy in the region a ^ r is presented in the form 



'a<r<b 



C{s)\s=-i = Ei% + El!'>+AE, 



(fe) 



(60) 



where E^"^^ (S^"^) is the vacuum energy for the region outside (inside) a cylindrical shell with 
radius a (b) and the interference term is given by the formula 



AE 



1^ oo 

(47r) 2 y-y / 



2 J n=0 
n Jn 



du (u 



m 



I 2 In 
ou 



llC\au)Kr{bu 



K';^\au)l)°> {bu 



ib). 



du u{u^ — rn^y 



In 



n=0 



lt\au)Kit\bu 



Kl^\au)l)::>{bu) 



ib), 



(61) 



To obtain the first of these formulae from the corresponding zeta function we have used the 
relation T{x) sinvrx = vr/r(l — x) for the gamma function. The interaction part of the vacuum 
energy (|61|) is negative for Dirichlet or Neumann boundary conditions and positive for Dirichlet 
boundary condition on one shell and Neumann boundary condition on the another. By the way 
similar to that used before for the case of the interaction forces it can be seen that in the limit 
a,b ^ oo for fixed b — a the corresponding result is obtained for parallel plates. In the limit 
a ^ for Aa,Ba 7^ the main contribution into the interaction part of the vacuum energy 
comes from n = term. By using the expansions for the Bessel modified functions for small 
values of the argument to the leading order one finds 



AE^ A 



D- 



TiAafl 
2Ba 



du u{u^ — m?) 2 



4\bu) ■ 



(62) 



In the same limit and for Dirichlet boundary condition on the inner cylinder, Ba = 0, the leading 
behavior for AE is obtained from (|62() by the replacement Aaa/Ba ln(a/6). In the case of 
Neumann boundary condition on r = a, the main contribution comes from n = and n = 1 
terms and AE vanishes as a^: 



^q2 roc 

AE^Ad— I duu^iu"^ 



m 



D-3 
I 2 



E' 



K^!'\bu) 



n=0 



r(b) 



(bu) 



In the limit 6 — > 00 and for a massless scalar field we have an asymptotic behavior with the 
leading term coming from n = summand: 



ttA 



D 



26^-1 ln(a/6) 



du u 



D-2 



Koju) 



(63) 



assuming that AajAf, ^ 0. For Aa = and Af^ ^ the main contribution into the interaction 
part of the vacuum energy comes from n = and n = 1 terms: 



na'^A 



2bD+ 



- 1 dun ^ - 

n=0 " 



Kn{bu) 



{bu) ■ 



(64) 



For the Neumann boundary condition on the outer cylinder, = 0, in the integrands of (|63() 
and (|64|) instead of ratio of the Bessel modified functions the ratio of their derivatives stands. 
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For a massive field and large values for the radius of the outer cylinder, under the condition 
mb S> 1 the main contribution into the integral over u in Eq. (|61() comes from the lower limit of 
the integral. By using the asymptotic formulae for the Bessel modified function for large values 
of the argument, to the leading order we find 



,1-D 



m 



' In. '(am) 



D-3 O-l 



(65) 



am] 



and the interaction part of the vacuum energy is exponentially suppressed. In figure we have 
plotted the dependence of the interaction parts (full curves) in the total vacuum energy on 
the ratio a/b for Z) = 3 massless scalar fields with Dirichlet, Neumann and Robin boundary 
conditions. For the Robin case we have chosen the parameters in the boundary conditions as 
Ba = 0, Bb/iAbb) = -0.2. 




2Ap(suif) / 



Figure 3: Interaction parts of the total and surface energies as functions of a/b. The graphs 
are plotted for D = 3 massless scalar fields with Dirichlet (D) and Neumann (D) boundary 
conditions and in the case of Robin (R) boundary condition with Ba = 0, Bh/{Af)b) = —0.2. 

To obtain the total Casimir energy, E, we need to add to the energy in the region between 
the shells, given by Eq. (|61|1 . the energies coming from the regions r ^ a and r ^ b. As a result 
one receives 

E= E^^^ + AE, (66) 

j=a,b 

where E^^^ is the Casimir energy for a single cylindrical shell with the radius j. The latter is 
investigated in Ref. j22j as a function on the ratio of coefficients in Robin boundary condition. 

(i) 

In particular, for Dirichlet and Neumann massless scalars one has the Casimir energies Ef) = 
6.15- 10~^j~^ and -E^^ = -1.42- 10~^j~^ (see ^H). Note that in the total vacuum energy single 
boundary parts dominate for small values of the ratio a/b and the interaction part is dominant 
for a/b < 1. In particular, combining the results for the single surface energies with the graphs 
from figure 13 we see that the total Casimir energy for a massless Dirichlet scalar in the geometry 
of two cylindrical shells is positive for small values a/b and is negative for a/b < 1. For Neuamnn 
case the vacuum energy is negative for all values a/b and has a maximum for some intermediate 
value of this ratio. 
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The total volume energy in the region a ^ r ^ 6 is derived by the integration of [j-component 
of the volume energy-momentum tensor over this region: 

l-b 

=2^ / drr(0|ro°|0). (67) 

J a 

Substituting mode-sum expansion into formula after the integration for the volume 
part of the vacuum energy we obtain 

i^^™') = i?a^.<6 - vr(4^ - 1) ^(0|'^'|0).=,. (68) 

j=a,b ^ 

As we see this energy differs from the total Casimir energy (|53|) (see Ref. |36j for the discussion 
in general case of bulk and boundary geometries). This difference is due to the presence of the 
surface energy located on the bounding surfaces. By using the standard variational procedure, 
in Ref. j36j it has been shown that the energy-momentum tensor for a scalar field on manifolds 
with boundaries in addition to the bulk part contains a contribution located on the boundary. 
For an arbitrary smooth boundary dMs with the inward-pointing unit normal vector n', the 
surface part of the energy-momentum tensor is given by the formula 



T^T"''^ =5{x-dMs)Tik (69) 



ik 

with 

Tik = i^^Kik - {2i - l/2)hik^n^Vi^, (70) 

and the 'one-sided' delta-function 5{x; dMg) locates this tensor on dMg. In Eq. ((7n|l. Kik is the 
extrinsic curvature tensor of the boundary dMg and /ij^ is the corresponding induced metric. In 
the region between the cylindrical surfaces for separate boundaries one has K^'j} = —jnjSfS'^, 
j = a, b. Now substituting the eigenfunctions into the corresponding mode-sum formula and 
using the boundary conditions, for the surface energy-momentum tensor on the boundary r = j 
one finds 



(0|r(^')'|0) = [-Un.dM + m - l/2)6fA,/B,\ (0|<^^|0),=„ (71) 
for z, fc = 0, 1, . . . , D, and {0\t[^^^\0) = 0. From ((ZU for the surface energy we obtain the formula 
= 27T ['drr{0\Tt''^^\0) = 7r{4^ - 1) V :^(0|(p2|o^ ^72) 

which, in accordance with (|68jl . exactly coincides with the difference between the total and 
volume energies. Of course, the VEVs of the field square on the right of this formula and, 
hence, the surface energy-momentum tensor diverge and need some regularization with further 
renormalization. Note that due to the surface divergences the subtraction of the boundary- 
free part is not sufficient to obtain the finite result and additional renormalization procedure is 
needed. In particular, the generalized zeta function method is in general very powerful to give 
physical meaning to the divergent quantities. However, in this paper we will not go into the 
details of the renormalization for the surface energy-momentum tensor for a single cylindrical 
boundary. This investigation will be presented in the forthcoming paper [301 ■ Here we note 
that after the subtraction of the boundary-free part the remained divergences in the surface 
energy-momentum tensor are the same as those for a single surfaces when the second surface is 
absent. The additional parts induced by the presence of the second surface are finite and can 
be obtained by using the representation of the VEV for the field square given by formula ((211) • 
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By taking into account the relation (|43|) . the surface energy on the boundary r = j is presented 
in the form 

^jsurf) ^^gurf)^^^j.urf)^ (73) 

where E^^^^^ is the surface energy for a single cylindrical boundary with radius j when the 
second boundary is absent and the term 

A R ^ , f'CXi D — 3 

AEf""''^ = -7t{4C-1)Ad^^Y1 / duu{u''-m'')~njr^{au,bu) (74) 



is induced by the presence of the second boundary. The latter is finite for all non-zero intersurface 
distances. Note that on the base of relation (|45|) it can also be written in the form 



n J m 

li''\au)Ki'\bu 



Bjn{u) dj 



n=0 

In 



(75) 



where Bjn{u) is defined by Eq. Now, by using formulae (|in)) . (|HT|). l(75|) it can be explicitly 

checked the relation 

= AS(™') + AEf'''^^ (76) 

j=a,b 

for the interaction parts of the separate energies. In figure 01 we have plotted the interaction 
part of the surface energy (dashed curve) for the case ol D = 3 massless Robin scalar with 
Ba = 0, Bf,/{Ahb) = —0.2 as a function on the ratio a/b. For the considered example this energy 
is located on the surface r = 6 — of the outer cylinder. The surface energies for Dirichlet and 
Neumann scalars vanish. 

Now let us explicitly check that for the interaction parts the standard energy balance equation 
is satisfied. We expect that in the presence of the surface energy this equation will be in the 
form 

^(surf) 

dE = -pdV + y dS(^), (77) 

^-^ 27ri 

j=a,b 

where V = ■K{b'^ - a^) and S^-') = 27rj are the volume and surface area per unit hypersurface in 
the axial direction. In Eq. H77jl . p is the perpendicular vacuum stress on the boundary and is 
determined by the vacuum expectation value of the ]^-component of the bulk energy-momentum 
tensor: p = — (0|Tj'^|0). From equation (|77|) one obtains 

^ = 2^jn,p(^) + -^^Vj, (78) 

with p^^^ being the perpendicular vacuum stress on the boundary r = j. Assuming that relation 
()78() is satisfied for the single boundary parts, for the interference parts we find 

^ = 2njn4l^+AEf-'^/j. (79) 

Now by taking into account expressions (|47|) . (|6H) . (|75j) for the separate terms in this formula 
and integrating by part in ()61() . we see that this relation indeed takes place. Hence, we have 
explicitly checked that the vacuum energies and effective pressures on the boundaries obey the 
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standard energy balance equation. Note that here the role of the surface energy is crucial and 
the vacuum forces acting on the boundary evaluated from the bulk stress tensor [determined by 
^(int)]' ™ general, can not be obtained by a simple differentiation of the total vacuum energy. 
The second term on the right of formula (|79j) corresponds to the additional pressure acting on 
the curved boundary. Noting that l^E^^^^^"^ /{2iTj) is the corresponding surface energy density, 
we see that this pressure is determined by Laplace formula. The total pressure on the boundary 
evaluated as the sum of bulk and surface parts is related to the total vacuum energy by standard 
formula and does not depend on the curvature coupling parameter. 

6 Conclusion 

In the present paper, we have investigated the one- loop quantum vacuum effects produced by two 
coaxial cylindrical shells in the + l)-dimensional Minkowski spacetime. The case of a massive 
scalar field with general curvature coupling parameter and satisfying the Robin boundary con- 
ditions on the boundaries is considered. To derive formulae for the VEVs of the field operator 
squared and the energy-momentum tensor, we first construct the positive frequency Wightman 
function. This function is also important in considerations of the response of a particle detector 
at a given state of motion through the vacuum under consideration 38 . The application of a 
variant of the generalized Abel-Plana formula to the mode-sum over zeros of the combinations 
of the cylindrical functions allowed us to extract the parts due to a single cylindrical boundary 
and to present the second boundary induced parts in terms of the exponentially convergent 
integrals. For the exterior and interior regions of a single cylindrical shell the Wightman func- 
tions are given by formulae (PUj) and pi]) respectively. The second boundary induced parts are 
presented by the last terms on the right of formulae (|19|) and (|23|) . The VEVs of the field square 
and the energy-momentum tensor are obtained by the evaluation of the Wightman function 
and the combinations of its derivatives in the coincidence limit of arguments. In both cases the 
expectation values are presented as the sum of single boundary induced and interference terms. 
The surface divergences in the VEVs of the local observables are contained in the single bound- 
ary parts and the interference parts are finite on both boundaries. In particular, the integrals 
in the corresponding formulae are exponentially convergent and they are useful for numerical 
evaluations. Due to the presence of boundaries the vacuum stresses in radial, azimuthal and 
axial directions are anisotropic. For the axial stress and the energy density we have standard 
relation for the unbounded vacuum. We have considered various limiting cases of the formulae 
for the interference parts. In particular, in the limit a, 6 — > oo for a fixed value b — a, we recover 
the result for the geometry of two parallel Robin plates on the Minkowski background. The 
vacuum forces acting on boundaries are considered in section 0J These forces contain two terms. 
The first ones are the forces acting on a single surface then the second boundary is absent. Due 
to the well-known surface divergences in the VEVs of the energy-momentum tensor these forces 
are infinite and need an additional regularization with further renormalization. The another 
terms in the vacuum forces are finite and are induced by the presence of the second boundary. 
They correspond to the interaction forces between the boundaries and are determined by for- 
mula (|44j) or equivalently by formula (|47j) . For the Dirichlet and Neumann scalars these forces 
are always attractive and they are repulsive for the mixed Dirichlet-Neumann case. In the case 
of general Robin scalar the interaction forces can be both attractive or repulsive in dependence 
of the coefficients in the boundary conditions and the distance between the boundaries. As an 
illustration, in figure |2l we present an example when the interaction forces are repulsive for small 
distances and are attractive for large distances. This provides a possibility for the stabilization 
of the radii by vacuum forces. However, it should be noted that for the reliable predictions 
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regarding quantum stabilization, the renormalized single shell parts also should be taken into 
account. In section [S] we consider the total vacuum energy in the region between the cylin- 
drical surfaces, evaluated as the sum of the zero-point energies for elementary oscillators. It 
is argued that this energy differs from the energy, obtained by the integration of the volume 
energy density over the region between the boundaries. We show that this difference is due to 
the presence of the surface energy located on the bounding surfaces. Further for the evaluation 
of the total and surface energies we use the zeta function technique. They are presented as the 
sum of single boundary and interaction parts. The latter are given by formula H61|) for the total 
vacuum energy and by formula (|75|) for the surface energy and are finite for all nonzero values 
of the intersurface separation. In the total vacuum energy single boundary parts dominate for 
small values of the ratio a/b and the interaction part is dominant for a/h < 1. For an arbitrary 
number of spatial dimensions and independent on the value of the mass, the interaction part 
of the vacuum energy is negative for Dirichlet or Neumann boundary conditions and is positive 
for Dirichlet boundary condition on one shell and Neumann boundary condition on the another. 
Further, we have shown that the induced vacuum densities and vacuum effective pressures on 
the cylindrical surfaces satisfy the energy balance equation (|77j) with the inclusion of the surface 
terms, which can also be written in the form H79|) . 
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